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Abstract 

We report on our first experiences with a mixed action setup with overlap valence quarks and 
non-perturbatively 0(a) improved Wilson sea quarks. For the latter we employ CLS Nf = 2 
configurations with light sea quark masses at small lattice spacings. Exact chiral symmetry allows 
to consider very light valence quarks and explore the matching to (partially quenched) Chiral 
Perturbation Theory (ChPT) in a mixed e/p-regime. We compute the topological susceptibility 
and the low-lying spectrum of the massless Neuberger-Dirac operator for three values of the sea 
quark mass, and compare the sea quark mass dependence to NLO ChPT in the mixed regime. This 
provides two different determinations of the chiral condensate, as well as information about some 
NLO low-energy couplings. Our results allow to test the consistency of the mixed-regime approach 
to ChPT, as well as of the mixed action framework. 

PACS numbers: 



I. INTRODUCTION 



Thanks to the theoretical and algorithmic improvements of recent years, and to the ever 
increasing computational power available, state-of-the-art Lattice QCD simulations now 
easily reach dynamical pion masses in the 200-300 MeV ballpark. 1 In this mass region the 
effective description of the dynamics of pseudo-Goldstone bosons at low energies by means 
of chiral perturbation theory (ChPT) is expected to work well at a quantitative level. This 
gives rise to a fertile interaction: by matching Lattice QCD and ChPT results it is possible, 
on the one hand, to test the effective description vs. the fundamental theory; and, on 
the other hand, low-energy constants (LECs) can be determined from first principles, thus 
providing a sounder foundation to phenomenological applications of ChPT. 

One particularly interesting aspect of the matching between QCD and ChPT is the role of 
finite volume effects [2H5] (we will always assume that the theory lives in an Euclidean four- 
volume V = L 3 x T). While for large enough values of L (one typical estimate is M n L > 4) 
the latter are expected to be strongly suppressed, and give rise only to small corrections to 
the infinite volume expansion in powers of pion momenta, the situation changes completely 
when the Compton wavelength of pions approaches L, i.e. M^L ~ 1. In this regime slow pion 
modes, strongly affected by the finite volume, dominate the path integral in the effective 
theory, and the expansion in powers of M^/A^ breaks down. A new power-counting for 
this e-regime was proposed in [3J, which implies a rearrangement of the chiral expansion, 
in which quark mass effects are suppressed relative to volume effects. As a consequence, 
less LECs appear at any given order in the expansion relative to the infinite volume case, 
and the approach to the chiral limit is therefore more universal. This regime leads to a 
very different setup for the determination of LECs, which offers both the potential to obtain 
cleaner computations of some of the latter (those whose effects are unsuppressed in the 
quark mass), and a cross-check of the systematic uncertainties of "infinite" volume studies. 

Another key property of finite volume chiral regimes (e-regime) is that the partition 
function of ChPT (at leading-order in the e-expansion) has been shown to be equivalent to a 
random matrix theory (RMT) jEH9], where many analytical predictions can be obtained for 
spectral quantities, such the spectral density or the distribution of individual eigenvalues. 

1 See e.g. the review p] at Lattice 2009. 
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These predictions are expected to be valid also for the spectrum of the Dirac operator in 
this regime, and have been tested both in quenched [T0hT2"] , Nf = 2 [TBTTTo] and Nf = 2 + 1 
[EH HZ] QCD. Since RMT predictions depend on one free parameter that corresponds to the 
chiral condensate (at the leading-order in the matching of ChPT and RMT), they provide 
yet another way of studying chiral symmetry breaking, using simple spectral observables. 

As we will see the matching of ChPT in the e-regime holds up to NLO in the chiral 
expansion, both for full and partially-quenched (PQ) situations as long as some quarks 
remain in the e-regime, because the path integral is dominated by the zero-modes of the 
lighter pions. The matching of QCD to RMT can therefore be extended to unphysical 
situations which are however more favorable from the computational point of view. 

Obviously enough, an adequate treatment of chiral symmetry on the lattice is especially 
relevant in this context. While simulations of Nf = 2(+l) QCD with full chiral symmetry 
have proven feasible, they are still limited to relatively small values of the inverse lattice 
spacing and/or physical volume [T5]. A way to overcome this is to use a mixed action 
approach [TjJJ [2U] , in which chiral symmetry is exactly preserved at the level of valence quarks 
only. Our aim is to develop such a framework by considering overlap valence quarks on top 
of Nf = 2 Coordinated Lattice Simulations (CLS) ensembles, 2 obtained from simulations 
with non-perturbatively 0(a) improved Wilson sea quarks. 

We will use this method to study the matching of QCD to ChPT and RMT in a mixed 
regime, in which sea quark masses are in the p-regime and valence quark masses are in the 
e-regime [2D 122]- From this matching we will be able to extract the Nf = 2 low-energy 
couplings £ and L 6 . 

Furthermore, the use of overlap fermions allows us to measure the topological suscepti- 
bility of the dynamical Wilson configurations. The dependence of this quantity on the sea 
quark mass has been derived in ChPT at NLO in infinite volume in [23] , and has been shown 
to depend on E and on a combination of L 6 , L 7 and L 8 . We will show how this prediction 
can be easily obtained in the mixed regime of ChPT and we will test it from the measured 
distribution of the topological charge. 

Obviously, mixed actions also have huge potential for phenomenological applications in 
which the exact preservation of chiral symmetry is greatly advantageous, e.g. to simplify 

2 https://twiki.cern.ch/twiki/bin/view/CLS/WebHome 
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the renormalisation of composite operators entering hadronic weak matrix elements. Along 
this line, first data for standard two- and three-point functions, as well as for correlation 
functions computed in the chiral limit via saturation with topological zero modes [21] 125] , 
will be covered in upcoming publications. 

The structure of the paper is as follows. In Section |TT] we review the main results from 
[2TJ [22] on the mixed-regime of ChPT and collect the results needed for our work, in which 
the sea quarks are degenerate and lie in the p- regime and valence quarks are in the e-regime. 
We will show that, up to NLO in the mixed-regime expansion, the partition functional 
matches onto a RMT, where the free parameter depends on the sea quark mass in a way 
that can be predicted from the matching. In Section |III| we present our numerical results for 
the low-lying spectrum of the overlap operator and the topological susceptibility on Wilson 



sea quarks and compare with the predictions of RMT and ChPT. In Section |IV| we present 
our results for the fits to ChPT predictions and extract the low-energy couplings. 



II. PROBING THE DEEP CHIRAL REGIME WITH MIXED ACTIONS 



A. Mixed chiral regimes 

While the exploration of chiral finite volume regimes ideally involves simulations with 
extremely light sea quarks, it is still possible to access them in a situation in which sea 
quarks have moderately larger masses. 

The first step is to formulate ChPT in a so-called mixed regime [21] [22], in which N s 
quarks have masses such that the p-regime requirement m s llV 3> 1 is satisfied, while N\ 
quarks have masses that fulfill the e-regime condition m\LV < 1. An appropriate power 
counting for this regime was first introduced in [2TJ [22] : 

Ptl ~0(e), L,T~0(l/e), ~ 0(e 4 ), m s ~0(e 2 ). (1) 

The partition function and meson correlators were computed to next-to-leading (NLO) order 
according to this power-counting, both in the context of fully dynamical quarks and also 
in various partially-quenched situations. The relevant PQ zero-mode integrals where also 
studied in [26]. In this work, we want to keep the sea quarks in the p-regime and the valence 
quarks in the e-regime. We briefly describe the results for this situation, and we refer to the 
original papers for details on the computations. 
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As is common in the e-regime of PQChPT, we need to consider sectors of fixed topological 
charge. We will be dealing therefore with the replica method in which the PQ limit is 
obtained in the limit Ni — > 0. 

The starting point is the parametrization of the Goldstone manifold according to the 
power-counting above. It can be shown that the non-perturbative zero-modes can be 
parametrized by a constant matrix, Uq G U(Ni) that together with various perturbative 
modes, £(x), T) ~ 0{e) span the full Goldstone manifold (£ contains all the non-zero momen- 
tum modes, while fj parametrizes the only perturbative zero-mode — for details see [2U 122])- 
The partition function in a sector of charge v is given by [2D 122] 

Z v ~ I [d£] [dfj] f [dU ] J(0 det(C7 )^exp ( - / d A xC(^U Q )\ . 

J JU{N t ) V J / 

(2) 

J(£) is the Jacobian associated with the parametrization [2TJ [221 121], an d both the La- 
grangian and the Jacobian can be perturbatively expanded in powers of £ and fj. The 
Lagrangian has an expansion in e of the form: 

£ = £ (4) + £ (6) + ..., (3) 

with terms up to C(e 4 ), up to £>(e 6 ), etc , while J(f) = 1 + C(e 2 ). 
The leading-order Lagrangian is found to be: 

2£ 



= Tr [d^d$ - |Tr [^MP^U, + U, 



Tr 

F 2 



V _ 



(4) 

where P; jS are the projectors on the light and heavy quark sectors and A4 is the quark mass 
matrix. At this order, there is a factorization of the perturbative and non-perturbative 
modes. 

An important observation from eq. Q is that the dependence on the non-perturbative 
modes, Uq, and on the light quark mass is identical to that of a theory with jVj quarks in the 
e-regime, but with the low-energy coupling, E, corresponding to a theory with Nf = Ni + N s 
flavours. This is to be expected since the heavier p-regime quarks behave, at the lowest 
energies, as decoupling particles that can be integrated out, but they are not quenched [22J. 

The quadratic form of the perturbative modes justifies their scaling with e. The last term 
in eq. (pj) could be treated as a perturbation as long as v ~ O(e ). However, in the partially- 
quenched case Ni = that we will be considering here, the distribution of topological charge 



is, on average, controlled by the sea quarks only. Indeed the v dependence of the leading- 
order function is found to be (after integrating over the perturbative modes) 



Z L V ° oc exp (-J^ J^p) J u{n) [dUo] det(U o yex V (|tt [l>lMI\(U + ul)fj 



(5) 



where 



Mi = (6) 

In the case Ni — > 0, the integral over the zero-modes is exactly one (as in the quenched case) 
therefore all the v dependence of the LO partition functional is in the explicit gaussian factor 
in eq. (j5]). The Leutwyler-Smilga result [2H] for the topological susceptibility is obtained: 

(v 2 ) = ^f 2 ^T ~ e~\ (7) 

a scaling that makes the last term in eq. Q of 0(e 4 ), and therefore of leading order. 

In fact we can easily push the computation of the topological susceptibility to NLO. In 
[29] [30] . the topological susceptibility has been computed in ChPT to NLO when all the 
quarks are in the p-regime. The same result should be obtained in the mixed-regime case 
when Ni — > 0, since only the dynamical p- regime quarks can contribute to the distribution of 
topological charges when Ni = 0. Indeed the NLO partition functional of ChPT in the mixed 
regime at fixed topology can be computed straightforwardly according to the power-counting 
rules of the mixed-regime, and the v dependence can be explicitly determined. It turns out 
that for Ni — > 0, all the ^-dependence comes from the integration of the perturbative modes 
£ and fj. The result for N s degenerate quarks is found to be 



2 m s El/ 
(v ) = 



1 N 2 s -l { M 2 SS , (M ss . , , , .,. 



2 



{Ll(p) + N.L r M + N s L r M) 



(8) 



where the function g%(M, L,T) contains the finite volume corrections to the closed pion 
propagator; its explicit definition can be found in [31]. This result agrees with the NLO 
results of (29] [30]. Note the appearance of L^(fi), for which no prediction has yet been 
obtained on the lattice. 
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In [6H9], it was shown that the e- regime zero-mode partition function at LO is that of a 
Random Matrix Theory (RMT) of matrices of size N, that depends only on the number of 
flavours, N, and the corresponding mass parameters rhi, viz. 



where Ai = d^rrii is the iVj x iVj mass matrix and there must be an identification Nrhi = 

From this relation, the microscopic spectral density of the Dirac operator, as well as 
higher order spectral correlation functions, can be related to those quantities computed in 
the corresponding RMT. Furthermore, the distribution of individual low-lying eigenvalues 
of the massless Dirac operator can also be predicted from this equivalence [32H31] , providing 
an efficient method to determine the chiral condensate, E. This relation has been tested in 
the quenched approximation and a good agreement has been found for volumes above 1.5 
fm or so [12]. In dynamical simulations, it has also been tested in [T3HT5] for Nf = 2 and in 
[T6| H7] for Nf = 2 + 1. More details on the RMT formulation will be given in the Section 



The rationale for expecting a matching of QCD to RMT relies on the existence of a 
regime where the chiral effective theory simplifies to a theory containing only the Goldstone 



in the mixed- regime above, there is a hierarchy of scales M vv <C L , which implies that 
we can integrate out the heavy scale L^ 1 to obtain a theory of zero-modes only, which we 
could call ZMChT (zero-mode chiral theory). We can obtain this theory from the full ChPT 
integrating the heavy modes order by order in the e-expansion. The difference between 
doing this matching in the e or the mixed-regime is the different assumption on the scaling 
of M ss . In the former case, M ss C L _1 and this scale is not integrated out (the ZMChT 
has therefore N + N s flavours), while in the latter L^ 1 ~ M ss and the zero-modes of the sea 
pions must be integrated out as well (the ZMChT has then N flavours). 

The matching of ChPT and ZMChT at LO in the mixed-regime can be easily derived from 
eq. Q: the ZMChT is simply ChPT at this order without the heavy modes (the integration 
over them gives an irrelevant normalization factor): 





he 



zero-modes, as depicted in Fig. [T] In fact, if we consider ChPT in the usual e- regime or 




(10) 
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FIG. 1: Chiral regimes of QCD showing the range of validity of the zero- mode chiral theory 
(ZMChT), which is equivalent to a RMT, and is obtained from ChPT by integrating out the non- 
zero momentum modes. A x ~ 4irF represents the chiral symmetry breaking scale. The LECs of 
the ZMChT (E e g) can be derived from the LECs of ChPT from matching, which can be done 
in two regimes depending on the size of the sea quark mass. For M ss <C 1/FL 2 , the e regime 
is appropriate, while for M ss > 1/FL 2 the mixed- regime has to be considered. The resulting 
ZMChTs have different number of flavours in the two cases. 



According to the eq. (J9J), this partition function is then equivalent to an Ni RMT. In partic- 
ular, it is important to stress that ZMChT has Ni flavours, while the full ChPT from which 
it is derived corresponds to Nf = Ni + N s flavours. In particular, for iVj — > 0, the ZMChT 
or RMT we expect to find is the quenched one, while the couplings should be those of an 
Nf = N s theory. 

The matching at NLO still does not modify the structure of the ZMChT theory. We 
could have anticipated this by realizing that there are no operators in the list of Gasser and 
Leutwyler that depend only on a constant Uq at 0(e 6 ). This does not mean however that 
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there are no corrections, simply that they can be absorbed in the couplings appearing at LO 



in eq. (10), that is E. Indeed at NLO, there are corrections to the zero-mode Lagrangian 



from the 0(p e ) terms 

£ (6) : 



. . . + — Tr 

y t 

16^Tt[P s M s P s }Tt 



PtMPilUo + U. 



(11) 



The integrations over the £, i] fields result in a change £ — > X e g- [2T 



yZMChTt 

£> v \NLO 



[dU }(det U yexp 



'U(Nt) 

with (for degenerate sea quarks) 



Tr 



Mi[U + W 



S c ff 1 _ 1 
~£~~ ~ F 2 

where 



16L 6 iV s M 2 - iV,A(0) - iV s A(M 2 s /2) + — A(0) 



N s 
NNi 



A(M 2 



M 2 = V ' y-' . 

v N f ss 



In dimensional regularization 
A(M 2 



M 2 / M 2 ' 

, -\ + ] og —)+g 1 (M,L,T), A(M 2 ) = A(M 2 



(4tt 



/// 



(12) 



(13) 



(14) 



(15) 



and A contains the expected UV divergence 

1 



A 



+ log4vr-7 + l-log/i 2 , e = 2-d/2, 



(16) 



that gets fully subtracted in the usual MS scheme. The small M expansion of A gives 



A(M 2 



1 M 2 
+ 



oo 1 

-A - 2lo gA iL) - £ r -^_ p n M^^L^\ 
' (n — 1 ! 



(17) 



KM 2 (4vr) 2 

v ' n=l 

where /3 n are the shape coefficients that depend only on the ratio T/L [3.1]. Note that the 
M — )■ limit of A(M 2 ) is well defined. 

In summary, up to NLO we have found that the ZMChT is equivalent to a RMT. Fur- 
thermore, the matching of this theory with ChPT gives the precise dependence on the sea 
quark mass of the coupling E e fj which is the only free parameter of the RMT theory. Testing 
this prediction will be one of the main results of this work. 

At this point it is interesting to discuss the possibility to have a smooth transition within 
the ZMChT regime between the Nf = Ni + N s effective theory and the iYj theory as the 



scale m s is increased. The authors of [H] have assumed that indeed this is possible and 
have found some evidence that the eigenvalue ratios seem to follow the dependence on m s 
predicted by the RMT or ZMChT. Such expectation would be justified if the conditions were 
such that FL 3> 1, because in this case the scale M ss can be neglected in the integration of 
the non-zero modes, as is done in the e-regime. However this is not true in practice, where 
FL ~ 1, and indeed even though the eigenvalue ratios in [H] showed roughly the dependence 
on m s predicted by the RMT (note that E e ff drops in the ratios), this is certainly not true for 
the eigenvalues themselves. The mixed-regime matching is the correct procedure to account 
for the correct m s dependence of E e g, for large enough m s . If FL is not sufficiently large, 
there is no warranty that the transition region (vertical band in Fig. [I]) can be modeled 
correctly by RMT. For a recent proposal to get predictions in the intermediate region from 
a resummation of zero-modes see 135 1. 



The distribution of the lowest lying eigenvalues of the Dirac operator is expected to match 
the prediction of RMT with Ni flavours and TV —> S cff V. S eff depends on the low-energy 
couplings of the Nf = Ni + N s theory. Note that the only NLO coupling entering is L 6 . 

Now since we want to consider the case of quenched light quarks, we have to take the 
limit Ni — > 0. S e ff has a finite replica limit given by 



M 2 

lim E e ff = £ { 1 + 



F 2 



02 logifjV^) N (M s , 

+ o„2*r + 1QN sL 6 (n) ~ j—^ log 



N s 8tt 2 N s b ^ J (4vr) 2 ^Vv^/i/ 



Pi N s 



N q F 2 VV 



-^ 9i (m ss /V2,L,t) I. (18) 



For the case T/L = 2, that we will be considering in our simulations, fi\ = 0.08360 and 
P 2 = —0.01295. For details on how to compute the shape coefficients we refer to Ref. [3T] . 
In Figure [2] we show the result of the ratio S e fj/E as a function of M 2 S /F 2 , for iV s = 2, 
F = 90 MeV and L r & (M p ) = 0.07 • 10~ 3 . The NLO corrections are quite significant, up to 
30% for the masses considered. 

Concerning the quenched limit of the zero-mode integral over U(Ni), a prescription using 
the supersymmetric or replica methods gives the well-known result for the quenched partition 
functional [361 137] that matches quenched RMT (qRMT). The low-lying eigenvalues of the 
Dirac operator should then follow the predictions of qRMT. Comparing the eigenvalues 
computed numerically in this PQ setup with the predictions of qRMT, we can extract S e ff 
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of eq. (18). If we do this for different values of the sea-quark mass, we can study the sea- 
quark mass dependence of E e g, from which we can in principle disentangle S and Lq(/j,), up 
to NNLO corrections (assuming we have an independent determination of F). 

A relevant question is however what are the eigenvalues that should be matched to RMT. 
Since there is a cutoff over which the ZMChPT should not be a good description, we expect 
that when the eigenvalues roughly reach such cutoff they should get significant corrections 
from the massive modes and therefore the matching to RMT should break up. A rough 
estimate would be the condition that A < mth, where m t h is the value of the quark mass 
corresponding to the p-regime. For example taking the value of m#, such that ML > 3, 
converts into the condition XLV < 9F 2 LT/2, which is roughly 6 — 7 for our lattices. This 
results in the expectation that only the few lowest eigenvalues (< 3 — 4 for our lattices) are 
below the threshold. For the largest eigenvalues , deviations from RMT could be sizeable, 
and the associated systematic uncertainty should be reduced by simulating at larger volumes. 

B. Random matrix theory 

We consider the gaussian chiral unitary model described by the partition function 

Ni 



Z v {m^ ...,m N 



■ f ) = t dWe-^ wfw ^Y[det(D + m), (19) 
J i=i 
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where 

/ w . 

D=\ , (20) 

and W is a complex rectangular matrix of dimensions (N + v) x iV. Here iV plays the role 
of the space-time volume, multiplied by a constant. We are interested in the large- N scaling 
limit at fixed /ij = rhiN. The partition function Z v then provides an equivalent description 
of the zero mode-chiral theory partition function in eq. Q PU with Ni flavours of mass 
rhi and fixed topological charge u, with the identification Nrhi = vn^LV = fii. 



If x is the k-th smallest eigenvalue of the matrix y/W^W, the probability distribution 
associated to the microscopic eigenvalue ( = Nx can be written as 

pK(C;W)= /Vi /W-- f rfa-i^(Ci,---,C^i,C;W), (21) 

with < d < • • • < Cfc-i < C The explicit form of w£(Ci> ■ ■ ■ > Cfc-ij C> {/■*}) * s known in the 
microscopic limit [221 US]. For instance, in the quenched case iVj = one has 

fc-i 

^(Ci,-..,a-i,a;{o}) = ^e-^ 4 a[lc M1 n s (Ci,G) 4 x (22) 

i=l fe— l>i>i>l 

^(^(Cfc, G)> s(C*> Ci)> • • • > s(Cfe, C*-i)> s (Cfc, Cfc-i), Cfe, • • • , Cfc), 



2(fc-l) 



with 



and 



*(Ci,0) = V^-CJ, (23) 
detA^ 



^2(^1, • • • , s n ) 



AH ' 



[AW]^^^), AW= J] ( 24 ) 

n>i>j>l 

where Jj are modified Bessel functions. The coefficient W£ can be fixed such that the prob- 
ability p v k is normalized to unity. There is an interesting property, called flavour-topology 
duality, which manifests itself at zero mass: p^{C', {0}) depends on the number of dynam- 
ical flavours and the topological charge only through the combination £ = Ni + \u\. The 
microscopic spectral density 

oo 

pS(C;{A*}) = X)rf(C;{A*}) (25) 



k=l 
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coincides by construction with the one computed in the ZMChT [38]. For instance, the 
quenched LO spectral density is given by 

P5 9 (C; o) = | [MO 2 - J» + x(C)J»-i(0] • (26) 

The equivalence can be extended to generic n-point density correlation functions. It is pos- 
sible to show that probability distributions of single eigenvalues can be defined also in the 
chiral effective theory by means of recursion relations involving all spectral correlators [39J . 
The clear advantage of RMT is that the probability distributions are computable in prac- 
tice, while in the chiral effective theory explicit expressions are missing. By assuming this 
equivalence holds for all spectral correlators, it is then legitimate to match the expectation 
values of the low eigenvalues of the massless QCD Dirac operator A& with the predictions of 
the corresponding RMT. 

We will be considering here a situation where two flavours of degenerate sea quarks have 
masses that are sufficiently large to be in the p-regime. In this case the sea quark mass does 
not appear explicitly in the ZMChT/RMT, as we have discussed. The latter corresponds to 
a theory with iVj — > light flavours, that is quenched RMT (qRMT). The sea quark mass 
dependence comes in only through E e g(M ss ) and can be predicted at NLO, as explained in 
the previous section. Therefore we expect 

(Cfc)qRMT = ^cs(M ss )\ Ni=0 V(X k )Q CD (M ss ), (27) 



where M ss is the sea pion mass, Y> e g(M ss )\ N[=Q is given in eq. (18), and expectations values 
are computed in RMT as 

(a>q R MT = / dCrf(C;0)C- (28) 

In this matching we assume that the QCD quark masses, the eigenvalues of the Dirac 
operator and the quark condensate are properly renormalised. The prediction for the ratio 
(Ck ) qRMT /id) qRMT * s parameter- free and can be compared directly with (Afc)Q CD /(A;)Q CD at 
any fixed M ss . 

On the other hand, if we consider ratios at different sea quark masses of the form 
(Afe)Q CD (Mi)/(Afe)Q CD (M 2 ), (Mi j2 are two different sea pion masses), we can assume that 
they can be matched to qRMT with appropriate values £ e fj(Mi i2 ) of the effective chiral 
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P = 5.3, c sw = 1.90952, V/a A = 48 x 24 3 


label 


K 


aM ss 


A^fg 


D 4 


0.13620 


0.1695(14) 


156 


D 5 


0.13625 


0.1499(15) 


169 


D 6 


0.13635 


0.1183(37) 


246 (D 6a : 159, D 6b : 87) 



TABLE I: Simulation parameters. 



condensate. Therefore 



(A fc )^ CD (M 1 ) (C k ) 



qRMT 



S eff (M 2 



S eff (M 2 



Ni=0 



Ni=0 



(29) 



It follows that information on the mass dependence of E e g, and hence on Lq, can be obtained 
from suitable eigenvalue ratios. 



III. RESULTS ON DIRAC SPECTRAL OBSERVABLES 

We have carried out our computations on CLS lattices of size 48 x 24 3 . The configurations 
have been generated with non-perturbatively O(a) improved Wilson fermions at /3 = 5.3 
and sea quark masses given by k = 0.13620,0.13625,0.13635 |4"0] . The simulations have 
been performed with the DD-HMC algorithm [H]; further details can be obtained in [42J. 
The lattice spacing has been determined to be a = 0.0784(10) fm in [40J, which implies 
that our lattices have physical size L ~ 1.88 fm and sea pion masses of 426, 377 and 297 
MeV, respectively. However, preliminary results from more precise determinations through 
different methods yield a ~ 0.070 fm jl3]. We will consider both values in our analysis. 

Following (30], we will refer to our three lattices as D4, D 5 and D 6 . It has to be noted that 
for the D 6 lattice we have two statistically independent ensembles, that we dub D 6a and D 6 t>- 
We have analyzed 246 Dq configurations, 169 D5 configurations and 156 D4 configurations; in 
all cases successive saved configurations are separated by 30 HMC trajectories of length r = 
0.5. In Table|T]we collect the simulation parameters. The sea pion masses in lattice units are 
taken from [12] for the lattices D4, D 5 , while for D 6 we performed an independent evaluation 
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FIG. 3: MC history and distribution of the topological charge in lattices D4 and D5. 



from the pseudoscalar correlator computed on 96 CLS configurations. The resulting value 
implies M^L = 2.84(9), which complies with the stability bound for the simulation algorithm 
derived in jH] . On the other hand, M^L < 3 implies sizeable finite volume effects in p-regime 
physics, which in the present work are accounted for within ChPT. 

On these configurations we have built the massless Neuberger-Dirac operator J45J H6] 

£n = ={l+7 5 sign(Q)}, (30) 

CL 

with 

Q = y s (aD w -l-s), a=-^—, (31) 

where _D\y is the Wilson Dirac operator. The parameter s governs the locality of and 
has been fixed to s = 0.4 for all our simulations. A discussion on the locality properties of 
the Neuberger-Dirac operator in our setup can be found in Appendix [Aj 

Our Neuberger fermion code is the same used in previous quenched studies p2J EHS)] , 
and is designed specifically to perform efficiently in the e- regime [51]. Our data analysis 
methods, including a discussion of autocorrelations in the observables under consideration, 
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FIG. 4: MC history and distribution of the topological charge in lattice Dg. 
are discussed in Appendix [Bj 



lattice 






D 4 


0.01(55) 


9.9(1.5) 


D 5 


-0.24(40) 


6.93(98) 


D 6 


0.62(24) 


3.36(47) 



TABLE II: Results for the expectation value of the topological charge and its square. 
A. Topological charge 

A first, immediate application of having constructed the Neuberger-Dirac operator 
on a given dynamical configuration is the determination of the topological charge of the 
latter by computing the index of Z?n, 

v = n + — n_ (32) 
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where n + (n_) is the number of zero modes of with positive (negative) chirality. 

In the upper part of Figs. [3]-[4]we show the Monte Carlo history of the topological charge 
for our three lattices. The topology sampling proceeds smoothly, although there are clear 
hints at the presence of sizeable autocorrelations (cf. Appendix [Bj . The histograms in the 
lower panels show the distribution of the measured topological charges, which qualitatively 
exhibits the expected Gaussian-like shape and width. This finding is consistent with the 
study reported in [52], since our computations take place at a value of the lattice spacing 
sufficiently larger than the threshold a < 0.05 fm below which topology is expected to exhibit 
freezing symptoms. 

In Table Owe quote our results for the expectation values (y) and (v 2 ). 



B. Low modes of the Dirac operator 

We have computed the 10 lowest eigenvalues of the Neuberger-Dirac operator on lattices 
D 4 , D 5 and D 6 by adopting the numerical techniques described in [5T] . 

The eigenvalues of appear in general in complex conjugated pairs and lie on a circle 
in the complex plane 

D N V> = 7^ 7==(l-e*). (33) 

CL 

In order to compare them with the predictions of RMT, we have computed the projection 

m 

A = = =^2(1 -cos0). (34) 

We have evaluated expectation values at fixed absolute value of the topological charge \u\. 



In Fig. 10 we show the bare eigenvalues for \u\ = 0, 1, 2. 



Since the matching with RMT involves the parameter E, it is useful to first consider ratios 



of eigenvalues. In our case, following eq. (27), the QCD ratios (\k) Ul / (\i) U2 can be directly 



matched with the qRMT predictions {Ck) Vl / '{(i) V2 ■ In Tables |Vrij |VIII| (App. [Bp we report 
the results for eigenvalue ratios involving the four lowest-lying eigenvalues and topologies 
\v\ = 0, 1, 2, together with qRMT predictions. It should be pointed out that the matching to 
RMT should work provided X^effV is not much larger that 1. For the lattice parameters 
we are considering, we set out cutoff at k < 4 for which the parameter is below 10. Since 
k = 4 is probably borderline, we will not include it in the extraction of S e jj though. 
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FIG. 5: Eigenvalue ratios at fixed \u\, normalized to qRMT predictions, for the indices (k,l) = 
(2,1), (3,1), (4,1), (3,2), (4,2), (4,3). 



The ratios at fixed topological charge are shown in Figs. [TTJ[T2j[T3](App.|B|; moreover, in 
Fig. [5] we report the ratios (^k) u / (^i) u normalized to the corresponding qRMT predictions, 
for \u\ = 0,1,2 and for several combinations k, I given in the bottom of the plot. This 
allows to appreciate clearly the precision and level of agreement with qRMT of each specific 
case. Finally, ratios at fixed k involving different topological sectors are presented in Fig. 



M](App.[B]). 

While the RMT prediction seems to work well for ratios not involving Ai, the ratios 
(Afc)/(Ai) exhibit somewhat more significant deviations. On the other hand, ratios between 
eigenvalues in different topological sectors follow well RMT predictions also in the case of Ai, 



as shown in Fig. [T4J albeit with larger errors. The data presented in this work do not allow 
for a full assessment of the systematics of these deviations, as this would require e.g. further 
values of the lattice spacing and/or physical volume. It is worth noting however that there 
are no clear differences between the three lattices, which we can take as an indication that 
corrections associated with relatively small values of m s are small. Concerning finite lattice 
spacing effects, having an estimate of the associated corrections in Wilson ChPT [20l [531 E3] 
would be welcome, although our value of the lattice spacing is quite small. Finally, one has 
to keep in mind that the impact of autocorrelations on statistical errors cannot be estimated 
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FIG. 6: The bare effective condensate £ e ff extracted from the matching in eq. (35 ), for k = 1, 2, 3, 4 



and \v\ = 0, 1, 2. The data for D4 have been shifted in the horizontal axis for better clarity. 

accurately for our ensembles. While we have attempted to stay on the safe side by quoting 
conservative errors that ought to include autocorrelations properly, it cannot be excluded 
that some errors are underestimated. Details are provided in App. [Bj 

C. Effective quark condensate 



In the spirit of the mixed regime ChiPT analysis, our data also allow to study the mass 



dependence of the effective condensate, cf. eq. (29). In Table IX we report the values of the 



bare effective condensate extracted from the matching 

(Cfc)qRMT 



Ks(M s 



V{X k ) 



V 

QCD 



(35) 



for k — 1, 2, 3, 4 and \u\ = 0, 1, 2. The results are shown in Fig. |6j where one can observe that, 
at fixed value of the sea quark mass, E e g does not depend on k and u within the statistical 
precision (with larger errors for k = 1). By averaging over k — 2, 3 and \v\ = 0,1,2 we 
obtain 



a 3 S cff = 0.00113(3)(4) (D 4 ), 
a 3 S cff = 0.00103(3)(4) (D 5 ), 
a 3 S cff = 0.00066(2)(5) (D 6 ). 



(36) 
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Lattice 


am 


om R (MS, 2 GeV) 


D 4 


0.00954(8) 


0.01366(23) 


D 5 


0.00761(7) 


0.01090(19) 


D 6 


0.00445(22) 


0.00637(33) 



TABLE III: Bare and renormalised (MS scheme at 2 GeV) PCAC sea quark masses. 



The first error is the statistical one, while the second uncertainty is a systematic effects 
estimated by adding the values for k — 1 in the average. We have checked that including 
k = 4 in the fit does not change the values within the statistical accuracy but decreases 
slightly the errors. 



Finally, in Fig. 15 (App. |BJ) we show the ratios defined in eq. (29) for k = 1,2,3, 4 and 
\v\ = 0, 1,2: they can be matched to the ratios of £ e fj at different quark masses. Here we 
can see that, as expected, those ratios do not depend on the topology and on k, within the 
statistical uncertainties. 



IV. FITS TO NLO CHIRAL PERTURBATION THEORY 

On the basis of the evidence presented in the previous section, now we assume that the 
matching to ChPT in the mixed-regime works in this range of sea quark masses and volumes, 
and try to extract the low-energy couplings from the sea-quark mass dependence of the two 
quantities S eff and (^ 2 ). 



The NLO predictions from ChPT are summarized in eqs. (18) and (|8j). As expected they 
depend on the two leading order LECs, S and F, but also on the 0(p 4 ) couplings L 6 , L7 
and L$. 

We first consider the topological charge distribution. The statistical error in this quantity 
is fairly large, but it is encouraging to see that there is a very clear dependence on the sea 
quark mass as shown in Fig. [7j We have fitted both to the full NLO formula in eq. (|8]), 
and to the linear LO behaviour. In either case m s is taken to be the PCAC Wilson mass 



renormalised in the MS scheme at 2 GeV, tabulated in Table III The results for the lattice D4 
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0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 

FIG. 7: (u 2 ) versus the sea quark mass omf s (2GeV). The smaller errors are the statistical ones 
and the largest include our estimate of autocorrelations. The dashed and solid lines correspond to 
the best LO and NLO ChPT fits respectively (for a = 0.078 fm). 



and D 5 are taken from [42], while we have computed that of D 6 , 3 using the renormalization 
constants and improvement coefficients of [55H59] . 

At LO the slope provides a direct measurement of S in the same scheme. At NLO we 
fit for £ and the combination [Lg + 2(Lg + Lj)] (M p ) after fixing \x = M p and rewriting 
M ss = 2m s 'E/F 2 . The value of F is fixed to 90 MeV; the systematic uncertainty related to 
this choice is estimated by varying F by ±10 MeV. 



In Table IV we show the results of the LO and NLO fits. In the case of the NLO, there is 
a slight difference when the scale is taken to be a = 0.078 fm (40] or the preliminary result 
a ~ 0.070 fm (43). We quote both. The x 2 is better for the NLO fit, but it is not possible 
to exclude the LO behaviour without decreasing our statistical errors. In physical units we 
get for a = 0.078 - 0.070 fm: 



S MS (2 GeV) = 262gj[Jj MeV - 287gJ{£{ MeV , (37) 



J (34)(5) 



3 



(35)(5) 



(36) (7) 



3 



where the first error is coming out of the fit and the second is the effect of changing F . There 
have been previous studies of the dependence on the topological susceptibility on the sea 
quark mass [23, E0H62J , but as far as we are aware the fits in these works have not included 
the NLO chiral corrections. 



We thank A. Jiittner for providing the necessary Wilson propagators for 
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Y, m (2 GeV)a 3 


[Ll + 2{L r ,+L r 7 )}{M p ) 


X 2 /dof 


LO 


0.00182(16) 




1.2 


NLO (a=0.078 fm) 


0.00112(48) 


0.0023(43) 


0.02 


NLO (a=0.070 fm) 


0.00106(44) 


0.0018(30) 


0.03 



TABLE IV: Results from the chiral fits to eq. (181) 



Let us now turn to £ e fr • In this case, the dependence on m s is expected starting at NLO 
in ChPT. The results in the previous section indicate that indeed there is a significant sea 
quark mass dependence in this quantity. We perform a two-parameter NLO fit, where we 
fix F and fit for £ and Lg(M p ). As a function of the MS (2 GeV) sea quark mass m s , we 
expect therefore: 



Zf s X eS (m s ) 



S 1 + 



2m s E 
Pi 

2F 2 V 1 / 2 



1 

2 16tt 2 



\og(M p V 1 ^) + 32L r 6 (M p 



1 



16n 



f m s T, 
2 ° g \F 2 M 2 




(38) 



where we need the scalar density renormalization factor in the MS scheme for the valence 
overlap fermions. We have obtained a rough estimate of this factor by matching our valence 
and sea sectors at a reference value of the pion mass, computed with mass-degenerate quarks, 
following the method of [63]. We have done this at the unitary point on lattice D5; choosing 
the latter instead of our lightest point D 6 allows to avoid sizeable finite volume effects in 
the determination of pion masses. 

The sea pion mass in lattice units is aM ss = 0.1499(15) [42], while for the mass of the 
valence pion at bare valence quark mass am = 0.020 we obtain aM vv = 0.153(5) with 63 D5 
configurations. 4 In order to obtain the renormalisation factor we then apply the matching 
condition 

1 overlap 



/MS 



rn 



m MS (2 GeV) 



M rcf 



Wilson 



(39) 



4 The relatively small error for this limited statistics is a result of the use of low-mode averaging [47] in the 
computation of the two-point function of the left-handed current, from which the mass is extracted. 
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FIG. 8: 5] e ga 3 as a function of the sea quark mass omf s (2GeV). The dashed and solid lines 
correspond to the best fit for £ and Lg(M p ) at F = 90 MeV, taking the scale to be 0.070 fm and 
0.078 fm respectively. 



at aM re f = aM ss . The renormalised MS mass of the sea quark mass for the D5 lattice can 



be read from Table |III[ and we obtain 

Zp(2 GeV) = 1.84(10), (40) 

where the error is dominated by the one in aM vv , i.e. in the determination of the unitary 
point, and can be much improved with a larger statistics in the valence sector. Obviously, 
several checks need to be done to ensure that this result is robust, such as checking the 
dependence on the reference pion mass, as well as on the sea quark mass. A careful study 
of renormalization is beyond the scope of this exploratory study, and will be performed in a 
forthcoming publication. 

With this estimate for the renormalisation factor and fixing F = 90 MeV, the result we 
obtain from the fit is, for a = 0.078 — 0.070 fm 

E m (2GeV)a 3 = 0.00102(18) - 0.00099(16) (41) 
Ll(M p ) = 0.0015(11) - 0.0010(7). (42) 

The quality of the fit is shown in Fig. [8j Although the fit is good, it would be desirable to 
have more sea quark masses and smaller statistical errors to assess the systematics of this 
chiral fit. Particularly useful would be to test the finite-size scaling. Translating to physical 
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units we have 



E MS (2GeV) = 255}X{ MeV - 280^$ MeV , (43) 



;(14)(1) 



'(16)(4) 



(14)(4) 



(16)(5) 



where the only systematic error that has been estimated is that associated to the change of 
F by ±10 MeV. 

This value of E is consistent with the one obtained from the topological susceptibility 
above, and both are in nice agreement with the alternative determination of [M], that 
extracted the condensate from the spectrum of the Wilson-Dirac operator on Nf = 2 CLS 
configurations at the same lattice spacing and sea quark masses, but in a larger physical 
volume. A number of recent determinations of E for Nf = 2 can be found in the literature 
(for a recent review see [65J). The matching to ChPT has been done in the p-rcgimc [66-68J, 
and also in the e-regime in [TH [T5J [TTJ EHHZE]- Our determination uses instead a PQ mixed 
regime and has been obtained with significantly finer lattices than the latter. Although our 
result lies in the same ballpark as many of these previous determinations, it is necessary to 
quantify the systematic uncertainties involved in our calculation. The dispersion of results 
for E existing in the literature is rather disturbing, and it is very important to do a proper 
job at evaluating the systematic uncertainties: finite a, systematics of the chiral fits and 
finite-size scaling. 



V. CONCLUSIONS 



We have implemented a mixed action approach to lattice QCD in which sea quarks 
are non-perturbatively 0(a) improved Wilson fermions, while valence quarks are overlap 
fermions. As a first application we have studied the spectrum of the Neuberger-Dirac oper- 
ator, as well as the topological susceptibility, in the background of dynamical configurations 
at a « 0.07 fm, for three values of the Wilson sea quark mass. The mixed-regime of ChPT 
provides predictions for these observables and their sea quark mass dependence in terms of 
various Nf = 2 chiral low-energy couplings: at the NLO, they are E, L 6 and the combina- 
tion Lg + 2(Lq + Lj). We find that these NLO predictions describe our data well, and the 
extracted LECs are in good agreement with previous determinations. 

This exploratory study obviously needs several important refinements to quantify sys- 
tematic errors in a reliable way. Different volumes should be considered to test the expected 
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finite-size scaling. Also, larger volumes will allow to augment the number of eigenvalues that 
can be safely matched to RMT predictions, which in turn will provide a definitive assessment 
of the associated systematic uncertainty. The lattice spacing is not known very precisely 
and an accurate determination is under way by various CLS groups. Obviously other (3 
values need to be considered to attempt a continuum extrapolation. Finally the effects of 
autocorrelations, that we have observed, would need larger statistics to ensure fully reliable 
statistical errors. 

Our results show that new PQ setups where sea and valence quarks may lie in different 
chiral regimes (p and e) are tractable (unphysical) regimes from which chiral physics can be 
extracted. Mixed actions are adequate to treat such regimes, and constitute an interesting 
approach for those applications where chiral symmetry plays an important role. 
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Appendix A: Locality properties of the Neuberger-Dirac operator 



The Neuberger-Dirac operator D N has been defined in eqs. (30), (31). While Q is 
ultralocal, sign(Q) — Q/ (Q^Q) 1 ^ 2 in general couples all space time points. As a consequence, 
the locality of the Neuberger-Dirac operator is not granted a priori. In [72] it was verified 
that locality is preserved in the quenched case, for values of the lattice spacing around and 
above the one we are considering now. We apply the method used there to our specific case. 

We analyze the effect of the sign operator on a localized field r\: 

ip(x) = sign(Q)7?(x) r} a (x) = 5 xy S al (Al) 

where y is some point on the lattice and a runs over the colour and the spin indices of the 
field. We evaluate the function: 

/(r) = max{||^(x)|| 2 | \\x - y||i = r} , (A2) 

where || . . . ||i denotes the so called "taxi-driver distance". It is clear that locality is re- 
covered in the continuum limit if /(r) decays exponentially, with rate proportional to the 
cutoff 1/a. To check that this indeed happens we fit the lattice data to Ae~ Br in a range 
fmin < t < r max . The upper limit has to be set because the inaccuracy with whom we cal- 
culate the overlap operator becomes bigger than the value of f(r) at large enough distances. 
The parameters of the simulation were set to calculate reliably HVK^)!! 2 at least down to 
values of about 1CT 16 (corresponding to r ~ 28). However, as can be seen in Fig. |9j our 
implementation of the overlap operator is more precise and the picture of a decaying expo- 
nential only breaks down at r ~ 38 where f(r) ~ 10 -22 . 

We have taken in all the cases r m j ra = 14 and r max = 28 also to compare with the results 



in [72]. The parameter s in eq. (31 ) can be varied to improve the locality properties of the 
Neuberger-Dirac operator. We collect the results of our fit for different values of s in Table 
|V} The quoted error is obtained through jackknife with bin size 1. Among the values we 
adopted in our test, the choice s = 0.4 yields the most satisfactory locality properties for 
this (3 value, as shown in lower panel of Fig. [9j This result is in agreement with the previous 
studies in the quenched case. We will therefore adopt s = 0.4 in our study. 
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s 


B 


10 4 A 


x 2 


0.2 


0.711 ±0.081 


0.3 ±0.1 


0.70 


0.4 


1.001 ±0.032 


3.4 ± 1.5 


0.45 


0.6 


0.975 ±0.031 


3.7 ± 1.2 


0.42 



TABLE V: Results of the fit for dynamical configurations after jackknife resampling. The results 
for S = 0.2,0.4,0.6 are based respectively on a set of 13, 21, and 20 configurations, respectively. 



Appendix B: Statistical error analysis 



Autocorrelations 



We have studied the presence of autocorrelations in our observables in various ways: 

• Integrated autocorrelation times have been estimated by using the methods described 
in [73] and [H]. In the first case, the summation window for the normalized autocor- 
relation function is fixed by setting the parameter S of [73J to S = 2, while in the 
second case we stop the summation when the normalized autocorrelation function is 
zero within one sigma. 

• The impact of changing the bin size in jackknife resampling of data has been assessed 
for all observables. 

• The impact of autocorrelations on statistical errors has been estimated directly with 
the techniques described in |73j for all observables (again with S = 2). 

In the case of the D$ lattice, we have studied autocorrelations in the Dg a ensemble only, as 
considering it together with the independent D 6 b ensemble would result in an underestima- 
tion of autocorrelation effects. 

Our primary observables are the topological charge v and expectation values of Dirac 
eigenvalues at fixed topology. In order to construct meaningful autocorrelation times for the 
latter, we consider their ratio with the RMT prediction for the expectation value of rescaled 
eigenvalues, i.e. (\k)u/ (C/c)BJvn> = (X'esV)^ 1 , which is ^-independent up to higher orders 



in the e chiral expansion. Results for autocorrelation times are given in Table |VI[ While 
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FIG. 9: Upper panel: Result of the fit of f(r) to Ae Br for s = 0.4. The blue line is the best fit 
to the data in the range 14 < r < 28. Lower panel: Comparison between various values of s. 

uncertainties on r int are remarkably large, due to the fact that measurements are performed 
only every 30 HMC trajectories, the values indicate that results coming from successive 
configurations are not completely decorrelated, especially in the case of the D 6 lattice. 

This is further reinforced by the analysis of the dependence of jackknife errors on the 
bin size. Again, this exercise is constrained by limited statistics, as the number of available 
configurations decreases rapidly with increasing \u\. Still, for topologies \u\ = 0,1,2 it is 
possible to have meaningful errors up to bin sizes of at least 5 configurations on lattices 
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\ 
I 


(v 2 ) 
\ u l 




'int 


'int 


' int ' int 


D 4 


12.0(6.2) 


8.0(5.2) 


1.47(58) 0.83(12) 


D 5 


6.4(3.4) 


3.6(1.7) 


1.7(0.7) 1.1(0.4) 


D 6 a 


4.9(2.5) 


3.8(1.8) 


2.9(1.3) 2.0(0.8) 





{^l)u/ (Cl)RMT> 

r (1) r. (2) 

' int ' int 


{M)v/ (C2)rMT> 

T (1) T (2) 
' int ' int 


{M)v/ (C3)rMT> 

r (1) r (2) 

' int ' int 


{M)v/ (C4)rmt> 
r {1) r {2) 

' int ' int 


D 4 
D 5 
D 6a 


0.404(94) < 0.5 
0.56(15) 0.89(27) 
2.7(1.3) 2.3(1.1) 


1.32(60) < 0.5 
1.45(63) 0.89(27) 
4.0(2.1) 3.2(1.5) 


1.56(77) < 0.5 
0.93(34) 0.89(27) 
7.0(3.9) 4.4(2.4) 


1.22(56) < 0.5 
1.22(54) 0.89(27) 
7.4(4.0) 5.2(3.2) 



TABLE VI: Estimated autocorrelation times for the topological charge and its square (upper panel), 
as well as for the ratios (Xk)u/ (Cfc)RMT> up to k = 4 (lower panel). Superindices (1) and (2) refer to 
computations following [73] and [H], respectively. T mt is provided in units of gauge configurations, 
i.e. it should be multiplied times 30 to convert to number of trajectories and times 15 to convert 
to MC time. 

D 5 and D 6 . On lattice D 5 , errors for (\k)\u\ and ratios (Xk)\ v \/ (\j)\ v \ exhibit little or no 
dependence on the bin size, with errors increasing by at most 20 to 30%. On the other 
hand, on lattice D 6 errors consistently increase with the bin size, and, in cases where there 
is enough statistics to avoid an early loss of signal, they tend to saturate around bin sizes 
of the order of 3-4, at which point they are between 30% and 70% larger than with bin 
size 1. Finally, the errors taking into account autocorrelations computed following [73] are 
consistent with jackknife errors with bin size 1 on D5, while on D$ they are systematically 
consistent with those around which the jackknife bin dependence stabilizes, and in some 
cases even slightly larger. 

Regarding the topological charge, where statistics allows to trace the bin size dependence 
up to much larger values, the increase in the error of (u) and (z/ 2 ) is much more marked 
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than for Dirac eigenvalues. Typically, jackknife errors stabilize for bin sizes between 5 and 
10 for all three lattices, at which point they are larger by as much as 60% with respect to 
those for bin size 1. The analysis a la [75] yields comparable errors. 

Our conclusion is that there is evidence that autocorrelations affect the topological charge 
in all three lattices, while spectral observables are affected by detectable autocorrelations 
on D$ only. As a simple recipe to stay on the safe side when including autocorrelation effect 
in the errors, we quote results for spectral observables from the analysis with jackknife bin 
size 1 for lattices D 4 and D 5 and bin size 3 for lattice D 6 . In the case of (y) and (z/ 2 ) we 
quote jackknife errors for the bin sizes at which they stabilize. It has to be stressed that we 
find no significative evidence that autocorrelations depend on the sea quark mass in some 
systematic way. 

2. Systematic errors in the computation of Dirac eigenvalues 

The numerical computation of eigenmodes and eigenvalues of the Hermitian non-negative 
operator -DJj-Dn has been performed with the techniques described in [oT] . The accuracy 
to which eigenmodes are computed is bound by an input parameter, which in our case has 
been set to 1% in lattices D 4 and D 5 and 5% in lattice D 6 . On the other hand, for each 
separate eigenmode an a posteriori estimate of the actual error in the computation of the 
eigenvalue is produced by the program. Usually, this estimate yields an error one order of 
magnitude smaller than the nominal accuracy parameter mentioned above. 

This systematic error should, in principle, be added in quadrature to the statistical error 
of (Afc) and observables derived thereof. We have estimated its impact, and found that it 
is completely negligible with respect to statistical errors in lattices D 4 and D 5 , both if it is 
computed using the estimates on each eigenvalue accuracy and in the much more pessimistic 
case in which a flat error associated to the 1% nominal precision is set. In the case of Dq, 
the uncertainty coming from numerical error estimates is again negligible, but a flat error 
set at 5% yields uncertainties comparable to the statistical ones. However, our experience 
shows that the estimates produced by the program are in the right ballpark, and conclude 
that setting a flat 5% uncertainty in D 6 observables would be a gross overestimate of the 
effect. Hence, we have opted for neglecting this source of error in our final results. 
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\v\ 


fc/i 


D 4 


D 5 


D 6 


qRMT 





2/1 


3.10(47) 


2.76(26) 


3.51(36) 


2.70 





3/1 


5.09(71) 


5.20(45) 


5.98(72) 


4.46 





4/1 


7.89(96) 


7.54(65) 


8.7(1.0) 


6.22 





3/2 


1.64(9) 


1.88(9) 


1.70(8) 


1.65 





4/2 


2.55(16) 


2.73(12) 


2.49(15) 


2.30 





4/3 


1.55(7) 


1.45(6) 


1.46(4) 


1.40 


1 


2/1 


1.98(10) 


2.08(13) 


2.26(10) 


2.02 


1 


3/1 


3.10(16) 


3.45(23) 


3.56(19) 


3.03 


1 


4/1 


4.30(23) 


4.80(33) 


4.85(30) 


4.04 


1 


3/2 


1.56(4) 


1.66(6) 


1.57(4) 


1.50 


1 


4/2 


2.17(6) 


2.31(7) 


2.15(8) 


2.00 


1 


4/3 


1.39(3) 


1.39(3) 


1.36(3) 


1.33 


2 


2/1 


1.98(15) 


2.00(10) 


2.04(8) 


1.76 


2 


3/1 


3.02(21) 


2.99(16) 


3.02(10) 


2.50 


2 


4/1 


3.99(28) 


3.93(20) 


3.94(14) 


3.23 


2 


3/2 


1.53(5) 


1.49(4) 


1.48(4) 


1.42 


2 


4/2 


2.02(7) 


1.96(6) 


1.93(7) 


1.83 


2 


4/3 


1.32(5) 


1.32(2) 


1.30(3) 


1.29 



TABLE VII: Results for eigenvalue ratios at fixed topology {Xk) u / {Xi) u compared with qRMT 
predictions (Cfc)7<Ci)"- 
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k 


M/V2I 


D 4 


D 5 


D 6 


qRMT 


1 


1/0 


2.20(31) 


1.87(24) 


2.12(29) 


1.75 


2 


1/0 


1.41(11) 


1.41(10) 


1.36(14) 


1.31 


3 


1/0 


1.34(7) 


1.24(8) 


1.26(10) 


1.19 


4 


1/0 


1.20(4) 


1.19(6) 


1.18(7) 


1.14 


1 


2/0 


2.55(41) 


2.69(32) 


2.76(35) 


2.45 


2 


2/0 


1.62(14) 


1.96(14) 


1.61(16) 


1.59 


3 


2/0 


1.51(11) 


1.55(10) 


1.40(11) 


1.37 


4 


2/0 


1.29(7) 


1.40(7) 


1.26(7) 


1.27 


1 


2/1 


1.16(13) 


1.44(14) 


1.30(11) 


1.40 


2 


2/1 


1.15(8) 


1.39(9) 


1.18(8) 


1.22 


3 


2/1 


1.13(8) 


1.25(7) 


1.12(6) 


1.15 


4 


2/1 


1.07(6) 


1.18(6) 


1.07(4) 


1.12 



TABLE VIII: Results for eigenvalue ratios (Xk) Ul / (^k) U2 compared with qRMT predictions 

<ar/<ar 2 - 
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k 


M 


a 3 S cff (D 4 ) 


a 3 X cff (D 5 ) 


a 3 S cff (D 6 ) 


1 





0.00135(18) 


0.00121(12) 


0.00089(11) 


2 





0.00118(8) 


0.00118(6) 


0.00068(6) 


3 





0.00119(5) 


0.00104(5) 


0.00066(5) 


4 





0.00107(3) 


0.00100(4) 


0.00063(3) 


1 


1 


0.00108(6) 


0.00113(9) 


0.00074(5) 


2 


1 


0.00110(4) 


0.00110(5) 


0.00066(3) 


3 


1 


0.00105(4) 


0.00099(4) 


0.00063(2) 


4 


1 


0.00101(3) 


0.00095(3) 


0.00061(1) 


1 


2 


0.00130(13) 


0.00110(7) 


0.00078(2) 


2 


2 


0.00116(7) 


0.00096(5) 


0.00068(3) 


3 


2 


0.00108(7) 


0.00092(4) 


0.00065(2) 


4 


2 


0.00105(5) 


0.00090(3) 


0.00064(2) 



TABLE IX: The bare effective condensate S e g extracted from the ratios (Cfc)qRMT/(^(-^)QCD) f° r 
k = 1,2,3,4 and \u\ = 0,1,2. 
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FIG. 10: The first 10 bare low eigenvalues of the massless Neuberger operator at fixed topology 
\v\ = 0, 1, 2, for the three lattices D4, D5, Y>§. 
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FIG. 11: Eigenvalue ratios {\k) u / {\i) v at fixed topological charge z/ = 0. The horizontal lines 
represent the qRMT prediction. 



34 



u 

5.5 

5 
4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 



jft/l" 



,3/1 



L 2/l 



i 1 ! 



4/2 



1 * i 



3/2 



1/2 



* I I 



4/3 



2/3 



1/3 



O.S 



3/4 



± i 5 



2/4 



5 * 5 



1/4 



I s 



FIG. 12: Eigenvalue ratios {Xk) u /(Xi) u at fixed topological charge |z^| = 1. The horizontal lines 
represent the qRMT prediction. 
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FIG. 13: Eigenvalue ratios {\k) v / {XiY at fixed topological charge |r/| = 2. The horizontal lines 
represent the qRMT prediction. 
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FIG. 14: Eigenvalue ratios at fixed k, for different topological sectors. The circles correspond to 
k = 1, the squares to fc = 2, the triangles to k = 3 and the diamonds to A; = 4. The horizontal 
lines represent the qRMT prediction. 
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FIG. 15: Ratios of eigenvalues at different sea quark masses, (\k) u (M\) / (\k) v {M2) , for k = 1, 2, 3, 4 
and =0, 1, 2. 
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